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The behavior of collective oscillations of a trapped boson-fermion mixed condensate is studied in 
the sum rule approach. Mixing angle of bosonic and fermionic multipole operators is introduced so 
that the mixing characters of the low-lying collective modes are studied as functions of the boson- 
fermion interaction strength. For an attractive boson-fermion interaction, the low-lying monopole 
mode becomes a coherent oscillation of bosons and fermions and shows a rapid decrease in the 
excitation energy towards the instability point of the ground state. In contrast, the low-lying 
quadrupole mode keeps a bosonic character over a wide range of the interaction strengths. For the 
dipole mode the boson-fermion in-phase oscillation remains to be the eigenmode under the external 
oscillator potential. For weak repulsive values of the boson-fermion interaction strengths we found 
that an average energy of the out-of-phase dipole mode stays lower than the in-phase oscillation. 
Physical origin of the behavior of the multipole modes against boson-fermion interaction strength 
is discussed in some detail. 
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Collective oscillation is one of the most prominent phenomena common to a variety of many-body systems. The 
CN| . realization of the Bose-Einstein condensates (BEC) for trapped Alkali atoms offers a possibility to study such 
phenomena of quantum systems under ideal conditions. Up to now the experimental as well as theoretical |3|-|5| 
studies of collective motions of BEC have been intensively performed. Quite recently the degenerate Fermi gas of 
trapped 40 K atoms has been realized 0, which motivates the study of collective motion also in Fermi gases ffl. These 
developments further encourage the study of possible boson-fermion mixed condensates of trapped atoms MM- Now 
the static properties stability conditions [jl3|,[l4]], and some dynamical properties 15 1^] of trapped boson- 
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fermion condensates have been investigated. In the present paper, we study the behavior of collective oscillations of a 
boson-fermion mixed condensate at T = for both repulsive and attractive boson-fermion interactions. We adopt the 
sum rule approach ||,^| that has proved to be successful in the studies of collective excitations of Bose condensates. 
For the mixed condensate, we introduce a mixing angle of the boson/fermion excitation operators so as to allow the 
^j" 1 , in- and out-of-phase oscillations of the bose and fermi condensates. 

In the sum rule approach we first calculate the energy weighted moments m p = — Eo) p \(j\F\Q)\ 2 of the 

relevant multipole operator F, where |j')'s represent the complete set of eigenstates of the Hamiltonian with energies 
Ej, and |0) denotes the ground state. The excitation energy is expressed as hu = (m^/mi) 1 / 2 which provides a 
useful expression for the average energy of the collective oscillation |||7]] . The moments are calculated from formulae 
mi = ±(0| [F\ [H,F]] |0) and m 3 = ±(0| [[F\ H], [H, [H,F]]] |0). We consider three types of multipole operators 
, which are defined by 

Ft = E fA) ± E f°w> ( a = M > D > Q) w 

i=l i=l 

where the functions f a are defined by /m(^) = r 2 for monopole, /o(r) = z for dipole, and /q(r) = 3z 2 — r 2 for 
quadrupole excitations. The indices b, f denote boson/fermion, Nb,Nf the numbers of bose/fermi particles, and 
± correspond to the in-phase and out-of-phase oscillation of the two types of particles. We actually take a linear 
combination of the form 

F a (r;9)=F+cos8 + F-sm9 (~ < 6 < (2) 

parametrized by the mixing angle 9. We study the value of 9 that minimizes the calculated frequency lo for each a. 

We consider the polarized boson-fermion mixed condensate in spherically symmetric harmonic oscillator potential. 
The system is described by the Hamiltonian 

N b ( -2 i i N » 1 N f ( rf 2 1 1 N " Nf 

ff=E | + + b E *<fr -**)}+ E f£ + \ + h T,T, - m ( 3 ) 

i=l I j=l j i=l L J i=l j=l 



1 



where we assume the same oscillator frequencies and masses for bosons and fermions for simplicity. The coupling 
constants g, h are the boson-boson/boson- fermion interaction strengths represented by the s-wave scattering lengths 
abf, and a^f as g — AitTi 2 a^b / m and h — 4Ttti 2 cibf /m. The fermion-fermion interaction has been neglected as the 
polarized system is considered. Following the standard calculation procedure the excitation frequencies are obtained 
as: 
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(ii) Quadrupole 



(iii) Dipole 



lE bb + lE bf + 2(E Mn + lE bb + i A') cos sin - A sin" 5 




gg(g) 

UJ 



E+ o + 2E^ o cos9sint 



(4) 



2E kin cos 9 sin t 



| A sin 2 
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E£ +2E^ cos9sint 



(5) 
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Here we dehned £± n = E b Mn ± E^, E± 



Q sin 2 9 



2N- cos 9 sin 



(6) 



±i?ho an<1 = Nb±Nf, where -E^f,/^ and ^ho' J/ are respectively 



the expectation values of the kinetic and oscillator potential energies for boson/fermion in the ground state. Boson- 
boson and boson- fermion interaction energies have been denoted by Ebb and The quantities A, A' and f2 are 
given in terms of the boson/fermion densities n>b(r), nf(r) in the ground state by 



d 3 rr 2 



dnt(r) drib(r) 



dr 
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dn b (r) dn f (r) 
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where £ = y/h/muJo- One may use the stationary condition of the ground state, 

2E+ n - 2E ho + 3E bb + 3E bf = 0, 2E^ a - 2E~ o + 3E bb + A' = (8) 
in order to eliminate in cq.(Q) the dependences on Ebb,Ebf and A'. The monopole frequency is then rewritten as 



uj m (9) 



E kin 



2E kin cos 9 sin ( 



2A sin 2 ( 



Ei 



2E ho cos 9 sin ( 



(9) 



We have checked that the Thomas-Fermi calculation of the ground state adopted below gives rise to a negligible 
difference if one evaluates either the expression (0) or (Q). 

We calculate the ground state energies and densities of the boson-fermion mixed system in the Thomas-Fermi 
approximation which is valid for gNb ^> 1 and Nf ^> 1 except around the surface region ]To| |l1^ , |l5| . We take 
harmonic oscillator length £ and energy Hloq as units, and define scaled dimensionless variables: the radial distance 
x = r/£, boson/fermion densities pbj{x) = nbj(r)£, 3 /Nbj, and chemical potentials fibj = 2p,bj /Tvu>o- We solve the 
coupled Thomas-Fermi equations, 



gN b pb{x) + x 2 + hNfPf(x) = /lb, [6ir 2 Nfp f (x)} 2/3 + x 2 + hN b p b {x) = /I) 



(10) 



where g = 2g/hiUo£ y 3 and h = 2h/fiLUo^ 3 . 

One of the most promising candidates for the realization of the mixed condensate is the potassium isotope system. 
Precise values of the scattering lengths are not well known at present and different values have been reported 17 -191. 
We take for the boson-boson interaction the parameters of the 41 K- 41 K system in [T^] and a trapping frequency of 
450Hz which gives g = 0.2. For the boson-fermion interaction we take several values in the range h/g = h/g = 
—2.37 ~ 3.2. It should be noted that the interaction strength can be controlled using Feshbach resonances [ po[ . 

We have performed a numerical calculation for Nb = Nf = 10 6 . In the ground state the fermions have a much 
broader distribution than bosons because of the Pauli principle. Fermions are further squeezed out of the center 
for a repulsive boson-fermion interaction [h > 0) . They will eventually form a shell-like distribution around the 



2 



surface of bosons for h/g > 1 and will be completely pushed away from the center (n/(0) = 0) at around h/g ~ 3. 
For an attractive boson-fermion interaction, on the other hand, the central densities of the bosons and fermions 
increase together. The system becomes unstable against collapse at around h/g = —2.37 due to the strong attractive 
boson-fermion interaction |Jl3 . 

Figure 1 shows the kinetic energy, the oscillator potential energy, and the interaction energy contributions to the 
ground state energy against the parameter h/g. The figure shows also the quantities A and which represent the 
contributions of the boson-fermion interaction to the multipole frequencies and (|^). One may notice that 

the fermionic kinetic- and potential-energy contributions are a few times larger than the bosonic contribution in the 
present system. It is noted that A takes large negative values at both large negative and positive regions of h/g. In 
the former region the bose and fcrmi density distributions become coherent due to the attractive interaction, and the 
radial integral in eq.(|7]) takes a large positive value. In the opposite case (h/g » 1), the same integral changes sign 
because the fermions are pushed away from the center by the repulsive boson-fermion interaction, thus giving a large 
negative contribution in the bosonic surface region. In the region < h/g < 1, on the other hand, the integral is 
slightly positive and A takes a small positive value. The quantity f2 follows the same trend as A, but the absolute 
value is much smaller than A, as the most important contribution to the integral comes from the surface region where 
r > £. 

Once the ground state parameters are obtained the frequencies (||)-(H) are minimized against 8. Usually, the sum 
rule approach predicts a strength-weighted average energy of eigenstates for a given multipolarity. The calculated 
energy coincides with the true excitation energy if the relevant strength is concentrated in a single state. By adopting 
the minimization procedure we simultaneously determine the character of the low-lying collective mode and the 
corresponding average energy. The character of the mode is given by the value of 9, for instance, 9 ~ 7r/4 for 
the bosonic- and — ir/4 for the fermionic-modes, 9 ~ for the in-phase oscillation and tt/2 for the out-of-phase 
oscillation of the two types of particles. As there are two kinds of particles involved in eq.(2), one would expect an 
emergence of two types of collective oscillations for each multipole. Another collective mode would have a character 
orthogonal to the low-lying mode as far as the phase relation of the two operators in (2) is concerned. In the present 
approach we calculate the frequency of the latter mode, the high-lying one, from eqs.(^)-(|^), by using the operator 
sin6>L — F~ cos9l for each a, where the mixing angle 6l is the one determined for the low- lying mode. 

Figure 2 shows frequencies of the lower (solid lines) and the higher (dashed lines) modes for (a) monopole, (b) 
quadrupole and (c) dipole cases as functions of h/g. The corresponding mixing angles 9 determined by the minimiza- 
tion procedure for the lower mode are plotted in Fig. 3 for each multipolarity as a function of h/g. Below we discuss 
the behavior of the frequencies cu a by defining three regions of h/g: (I) h/g < 0, (II) < h/g < 1, (III) 1 < h/g. 

ajmonopole: 

For a non-interacting boson-fermion system the low-lying monopole mode is the fermionic oscillation with frequency 
uj m = 2wo, while the higher mode is the bosonic one with uj m = V§uo in the Thomas-Fermi approximation. Around 

h ~ one may obtain uj m ~ 2w (l + AhN* ) , uj m ~ VEluq (l - A'hg~iN^, where A = 3 • 6^ /4\/27r 2 and A' = 

7\/3/207r 2 -(87r/15)i The mixing angle for the lower mode is given as 8 M = -tt/4-<5m with 5 M = (5-65/V27r 2 )WVi 
This behavior is seen in region (II) where the boson-fermion interaction is weakly repulsive and the lower monopole 
mode is of a fermionic character, see Fig. 3 (solid line). Bosonic oscillation in this region is more rigid than the 
fermionic one because of the repulsive interaction among bosons. In region (I), the situation is quite different: The 
low-lying mode becomes a coherent boson-fermion oscillation as represented by the large negative value of A, and 
the excitation energy shows a sharp decrease towards the instability point h/g ~ —2.37 of the ground state (li), 
although ljm does not become exactly zero within our approximation. In this region the attractive boson-fermion 
interaction is much more effective in the excited state than in the ground state and cancels out the increase in the 
kinetic energy. In region (III), too, we find that the low- lying mode becomes an in-phase oscillation. Here, the boson 
and the fermion densities in the ground state are somewhat separated, and the in-phase motion which keeps this 
separation is energetically more favorable than the out-of-phase motion as seen in the value of A. 

b) quadrupole: 

For the quadrupole excitation (Fig. 2(b)), the lower (higher) energy mode is almost the pure bosonic (fermionic) 
oscillation over the broad range of the h/g values studied, Fig. 2 (dashed line). To the first order in h the frequencies 

of the lower and the higher quadrupole modes are given by ujL ~ V2cj (l-BhNe^, a;** ~ 2w (l - B'hgiN™), 

where B = 6*/4\/2V and B 1 = (15/8tt)s/(14- 6iV2w 2 ). The corresponding mixing angle for the lower mode is given 
by ®Q = 7r /4 + 5q with Sq = 2 3 B'hg~ N^o . For the quadrupole mode similar mechanisms as for the monopole mode 
are at work concerning the dependence on \h/g\. The role of the boson-fermion interaction is however much reduced 
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compared with the monopole case as seen by the factor 2/5 of eq.(^), which reflects that the quadrupole oscillation 
has five different components. Thus the quadrupole mode obtains an in-phase character only at large values of \h/g\. 
In the other region of \h/g\ the low- lying mode becomes a simple bosonic oscillation. This is because the fermionic 
mode costs a larger kinetic energy and favors 9 — ir/4 as seen in the term (E^ in — E[ in ) sin cos 9 in eq.(||). 

c)dipole: 

General arguments [Q show that for a harmonic oscillator external potential a uniform shift of the ground state 
density generates an eigenstate of the system, corresponding to the boson-fermion in-phase dipole oscillation with 
frequency Uq, This is evident in Fig. 2c) and also in eq.([j|) at 6 = 0. In the regions (I) and (III) the out-of-phase 
oscillation is unfavorable due to the same reason as for the monopole mode: It loses the energy gained in the ground 
state boson-fermion configuration. For a weakly repulsive h an interesting possibility arises: In the region (II) the 
out-of-phase mode of the boson-fermion oscillation lies lower than the in-phase mode. Let us first note that at h = 
the out-of-phase oscillation frequency becomes degenerate as the in-phase one because the bosonic and the fermionic 
dipole modes are independent. One may note that at h/g ~ I again the degeneracy occurs. Here the potential term for 



the fermion becomes almost linear to the fermion density itself (see, eq.(IO)). This suggests that the fermion density 
is determined almost entirely by the chemical potentials and becomes nearly constant as far as the boson density is 
finite. A uniform dipole shift of fermions thus produces almost no effect on bosons, and results in the degeneracy of 
the frequency. Between h/g — and I, the boson-fermion repulsion is weaker for the out-of-phase oscillation than the 
in-phase one (and hence the ground state) as reflected in the sign of f2. 

In the present paper, we studied collective oscillations of trapped boson-fermion mixed condensates using sum rule 
approach. We introduced a mixing angle of bosonic and fermionic multipole operators so as to study if the in- or 
out-of-phase motion of those particles is favored as a function of the boson-fermion interaction strength. For the 
monopole and quadrupole cases, the coupling of the bose- and fermi-type oscillation is not large for moderate values 
of the coupling strength h. At large values of \h/g\, the low-lying modes become an in-phase oscillation of bosons and 
fermions. This is especially so for the monopole oscillation at attractive boson-fermion interaction: The excitation 
energy of this mode almost vanishes at the instability point of the ground state. In the case of the dipole mode, in 
contrast, the in-phase oscillation remains an exact eigenmode with a fixed energy for harmonic oscillator potentials, 
while the average energy of the out-of-phase oscillation is strongly dependent on the boson-fermion interaction. Wc 
found that at weak repulsive values of the interaction the out-of-phase motion stays lower than the in-phase oscillation. 
In this paper we calculated also the frequencies of the high-lying modes for each multipole, by adopting the operators 
orthogonal to the low-lying modes. These modes, too, are collective in character and, in the present framework, their 
average frequencies showed rather strong dependences on the boson-fermion coupling strength. Deeper insight into the 
collective modes studied in this paper will require a detailed investigation of the solutions of, e.g., the self-consistent 
RPA type equations that allow an arbitrary radial dependence of the excitation operators. Studies in this direction 
are now in progress. 
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FIG. 1. (a) Ground state expectation values of the fermion kinetic energy E[ in (long dashed line), oscillator potential energies 
-E^o (short-dash-dotted line) and E^ o (long-dash-dotted line), boson-boson and boson- fermion interaction energies Ebb (dashed 
line) and Ebf (dotted line), and the quantity A (solid line) against the interaction strength ratio h/g. The values are given in 
the unit of Nhuio and are dimensionless. (b) The quantity Q in the same unit. 

FIG. 2. Excitation frequencies of collective oscillations ( a) monopole, b) quadrupole, c) dipole ) as functions of h/g — h/g, 
calculated based on the eqs.(Q)-(^|). The mixing angle 9 has been determined so as to minimize to for each multipole operator. 
The ordinates are given in the unit of wo and are dimensionless. The solid (dashed) lines are the lower (higher) energy mode 
for each oscillation. 

FIG. 3. Mixing angles of in/out-of-phase excitation modes which are determined by minimizing excitation energies in Fig.3. 
The solid line shows the monopole, the dashed line the quadrupole, and the dotted line the dipole oscillations. Shaded areas 
show mainly in-phase (hatches) and mainly out-of-phase (cross-hatches) regions in 8 Angles for pure bosonic- or fermionic-modes 
are also indicated. 
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FIG.l Ground state expectation values 




FIG. 2 Excitation frequencies 



